Physics:Electric-field integral equation

The electric-field integral equation is a relationship that allows the calculation of an electric field (E) generated by an electric current distribution (J).

Derivation

When all quantities in the frequency domain are considered, a time-dependency [math]\displaystyle{ e^{jwt} }[/math] that is suppressed throughout is assumed.

Beginning with the Maxwell equations relating the electric and magnetic field, and assuming a linear, homogeneous media with permeability [math]\displaystyle{ \mu }[/math] and permittivity [math]\displaystyle{ \varepsilon\, }[/math]: [math]\displaystyle{ \begin{align} \nabla \times \mathbf{E} &= -j \omega \mu \mathbf{H} \\[1ex] \nabla \times \mathbf{H} &= j \omega \varepsilon \mathbf{E} + \mathbf{J} \end{align} }[/math]

Following the third equation involving the divergence of H [math]\displaystyle{ \nabla \cdot \mathbf{H} = 0\, }[/math] by vector calculus we can write any divergenceless vector as the curl of another vector, hence [math]\displaystyle{ \nabla \times \mathbf{A} = \mathbf{H} }[/math] where A is called the magnetic vector potential. Substituting this into the above we get [math]\displaystyle{ \nabla \times (\mathbf{E} + j \omega \mu \mathbf{A}) = 0 }[/math] and any curl-free vector can be written as the gradient of a scalar, hence [math]\displaystyle{ \mathbf{E} + j \omega \mu \mathbf{A} = - \nabla \Phi }[/math] where [math]\displaystyle{ \Phi }[/math] is the electric scalar potential. These relationships now allow us to write [math]\displaystyle{ \nabla \times \nabla \times \mathbf{A} - k^2\mathbf{A} = \mathbf{J} - j \omega \varepsilon \nabla \Phi }[/math] where [math]\displaystyle{ k = \omega \sqrt{\mu \varepsilon} }[/math], which can be rewritten by vector identity as [math]\displaystyle{ \nabla (\nabla \cdot \mathbf{A}) - \nabla^2 \mathbf{A} - k^2 \mathbf{A} = \mathbf{J} - j \omega \varepsilon \nabla \Phi }[/math]

As we have only specified the curl of A, we are free to define the divergence, and choose the following: [math]\displaystyle{ \nabla \cdot \mathbf{A} = - j \omega \varepsilon \Phi \, }[/math] which is called the Lorenz gauge condition. The previous expression for A now reduces to [math]\displaystyle{ \nabla^2 \mathbf{A} + k^2\mathbf{A} = -\mathbf{J}\, }[/math] which is the vector Helmholtz equation. The solution of this equation for A is [math]\displaystyle{ \mathbf{A}(\mathbf{r}) = \frac{1}{4 \pi} \int \mathbf{J}(\mathbf{r}^{\prime}) \ G(\mathbf{r}, \mathbf{r}^{\prime}) \, d\mathbf{r}^{\prime} }[/math] where [math]\displaystyle{ G(\mathbf{r}, \mathbf{r}^{\prime}) }[/math] is the three-dimensional homogeneous Green's function given by [math]\displaystyle{ G(\mathbf{r}, \mathbf{r}^{\prime}) = \frac{e^{-j k \left|\mathbf{r} - \mathbf{r}^{\prime}\right|}}{\left|\mathbf{r} - \mathbf{r}^{\prime}\right|} }[/math]

We can now write what is called the electric field integral equation (EFIE), relating the electric field E to the vector potential A [math]\displaystyle{ \mathbf{E} = -j \omega \mu \mathbf{A} + \frac{1}{j \omega \varepsilon} \nabla (\nabla \cdot \mathbf{A})\, }[/math]

We can further represent the EFIE in the dyadic form as [math]\displaystyle{ \mathbf{E} = -j \omega \mu \int_V d \mathbf{r}^{\prime} \mathbf{G}(\mathbf{r}, \mathbf{r}^{\prime}) \cdot \mathbf{J}(\mathbf{r}^{\prime}) \, }[/math] where [math]\displaystyle{ \mathbf{G}(\mathbf{r}, \mathbf{r}^{\prime})\, }[/math] here is the dyadic homogeneous Green's Function given by [math]\displaystyle{ \mathbf{G}(\mathbf{r}, \mathbf{r}^{\prime}) = \frac{1}{4 \pi} \left[ \mathbf{I}+\frac{\nabla \nabla}{k^2} \right] G(\mathbf{r}, \mathbf{r}^{\prime}) }[/math]

Interpretation

The EFIE describes a radiated field E given a set of sources J, and as such it is the fundamental equation used in antenna analysis and design. It is a very general relationship that can be used to compute the radiated field of any sort of antenna once the current distribution on it is known. The most important aspect of the EFIE is that it allows us to solve the radiation/scattering problem in an unbounded region, or one whose boundary is located at infinity. For closed surfaces, it is possible to use the Magnetic Field Integral Equation or the Combined Field Integral Equation, both of which result in a set of equations with improved condition number compared to the EFIE. However, the MFIE and CFIE can still contain resonances.

In scattering problems, it is desirable to determine an unknown scattered field [math]\displaystyle{ E_{s} }[/math] that is due to a known incident field [math]\displaystyle{ E_{i} }[/math]. Unfortunately, the EFIE relates the scattered field to J, not the incident field, so we do not know what J is. This sort of problem can be solved by imposing the boundary conditions on the incident and scattered field, allowing one to write the EFIE in terms of [math]\displaystyle{ E_{i} }[/math] and J alone. Once this has been done, the integral equation can then be solved by a numerical technique appropriate to integral equations such as the method of moments.

Notes

By the Helmholtz theorem a vector field is described completely by its divergence and curl. As the divergence was not defined, we are justified by choosing the Lorenz Gauge condition above provided that we consistently use this definition of the divergence of A in all subsequent analysis. However, other choices for [math]\displaystyle{ \nabla\cdot\mathbf{A} }[/math] are just as valid and lead to other equations, which all describe the same phenomena, and the solutions of the equations for any choice of [math]\displaystyle{ \nabla\cdot\mathbf{A} }[/math] lead to the same electromagnetic fields, and the same physical predictions about the fields and charges are accelerated by them.

It is natural to think that if a quantity exhibits this degree of freedom in its choice, then it should not be interpreted as a real physical quantity. After all, if we can freely choose [math]\displaystyle{ \nabla\cdot\mathbf{A} }[/math] to be anything, then [math]\displaystyle{ \mathbf{A} }[/math] is not unique. One may ask: what is the "true" value of [math]\displaystyle{ \mathbf{A} }[/math] measured in an experiment? If [math]\displaystyle{ \mathbf{A} }[/math] is not unique, then the only logical answer must be that we can never measure the value of [math]\displaystyle{ \mathbf{A} }[/math]. On this basis, it is often stated that it is not a real physical quantity and it is believed that the fields [math]\displaystyle{ \mathbf{E} }[/math] and [math]\displaystyle{ \mathbf{B} }[/math] are the true physical quantities.

However, there is at least one experiment in which value of the [math]\displaystyle{ \mathbf{E} }[/math] and [math]\displaystyle{ \mathbf{B} }[/math] are both zero at the location of a charged particle, but it is nevertheless affected by the presence of a local magnetic vector potential; see the Aharonov–Bohm effect for details. Nevertheless, even in the Aharonov–Bohm experiment, the divergence [math]\displaystyle{ \mathbf{A} }[/math] never enters the calculations; only [math]\displaystyle{ \nabla\times\mathbf{A} }[/math] along the path of the particle determines the measurable effect.

References

• Gibson, Walton C. The Method of Moments in Electromagnetics. Chapman & Hall/CRC, 2008. ISBN:978-1-4200-6145-1
• Harrington, Roger F. Time-Harmonic Electromagnetic Fields. McGraw-Hill, Inc., 1961. ISBN:0-07-026745-6.
• Balanis, Constantine A. Advanced Engineering Electromagnetics. Wiley, 1989. ISBN:0-471-62194-3.
• Chew, Weng C. Waves and Fields in Inhomogeneous Media. IEEE Press, 1995. ISBN:0-7803-4749-8.
• Rao, Wilton, Glisson. Electromagnetic Scattering by Surfaces of Arbitrary Shape. IEEE Transactions on Antennas and Propagation, vol, AP-30, No. 3, May 1982. doi:10.1109/TAP.1982.1142818